Initial stages of thin film growth in the presence of island-edge barriers 
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qq , A model of submonolayer thin film growth is studied, where the attachment of atoms to island 

' edges is hindered by an energy barrier. A novel behavior of the density of islands, N 3 , is predicted 

^\ | as a function of flux F and temperature T. For example, N s scales as F x with X = 2i* /(i* + 3), 

where i* is the critical island size, in contrast with the standard result X = i* /{i* + 2). The theory 
' is applicable to surfactant mediated growth and chemical vapor deposition. It explains recent 

Oh. experiments, which are inconsistent with the standard theory. 
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Thin film growth processes are of tremendous importance for the fabrication of nanostructures and electronic 
devices. Technological advances in device miniaturization depend largely on the degree to which one can control 
the growth process, the epitaxial quality of the film and the morphology of the surface. It is therefore essential to 
understand the microscopic processes involved in thin film growth and especially their effect on the structure of the 
1 1 , film. 

Of particular interest are the initial stages of growth or the submonolayer regime, which is relatively easy to inves- 
tigate both experimentally and theoretically. It is possible to learn from such studies about the relevant microscopic 
processes and their respective energy barriers. For example, when the film evolves by nucleation, growth and coales- 
cence of two dimensional (2D) islands, the shape and size distributions of the islands and their dependence on physical 
parameters such as temperature, flux and coverage, yield information about various diffusion processes, attachment 
and detachment of atoms to and from island edges, etc. |Q. 

The first theories of diffusion of atoms in the presence of steps assumed that step edges are perfect sinks for 
adatoms Q . This assumption was later relaxed || , and kinetic coefficients were introduced to take into account the 
finite energy barriers || associated with attachment and detachment of adatoms to and from the edges. The perfect 
l— ~~ ', sink assumption has been shown to yield reasonable results under many experimental conditions. It was realized, 
however, that it fails in several important cases. For example, a significant barrier for attachment of adatoms to 
, steps from the terrace above, leads to a kinetic instability of the flat surface and to the growth of large mounds 
in ' H- An asymmetry in the barriers for adatom attachment from below and above the step can also lead to peculiar 
step bunching j(| as well as fingering [Q instabilities. Island-edge barriers may be important in surfactant mediated 
growth, where a surfactant atom can bind to an island edge. In order for an adatom to attach to the edge, surfactant 
atoms have to be removed from there. Kandel and Kaxiras explained ^,|) experimental results related to surfactant 
' mediated growth by assuming that the energy barrier for such a change in the atomic configuration is large [fio|| . 
^\ , Island-edge barriers may also occur in chemical vapor deposition (CVD). For example, during CVD of Si on Si using 
disilane, the surface is covered with hydrogen, which can bind to island edges similarly to surfactant atoms 

In this work, the effect of island-edge barriers on submonolayer growth is studied in the framework of rate equation 
theory [p], p^|Jl3| . In particular, the density of 2D islands on the surface, N S (F, T, 6), is calculated as a function of flux 
, F, temperature T and coverage 9. It is shown that island-edge barriers have a dramatic, experimentally observable, 
effect on the behavior of N s . Although finite island-edge barriers have been studied before @,|5| in the context of 
submonolayer growth, it is the first time their detailed effect on the dependence of the island density on flux and 
O ■ temperature is addressed. 

^ , The simplest scheme for the calculation of N s is the critical island approximation, where it is assumed that islands 
that contain more than i* atoms are stable, while smaller ones are not and can decay jlfj]. Within this scheme, N s is 
the density of stable islands, and a detailed balance relation is assumed to hold [jl7| between the densities of unstable 
' islands, N{ (i < i*), and the average adatom density n: 

QNi = (Qnye 0Ei , (1) 

where f2 is the atomic area of the solid, (3 = l/fc^T and Ei is the binding energy of an island of i atoms. 
Now one can write down the rate equation for the density of stable islands: 
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where C and I are the coalescence and nucleation rates per unit area. When the coverage is small, coalescence does 
not occur and C«0. Since this is the limit of interest in this work, C is completely neglected below. The nucleation 
rate is 

I = anN*S* . (3) 
Here a is the capture coefficient of a critical island, N* = Ni* and 

S* = Me-^ +E ^ . (4) 

In the equation above, Ed is the diffusion barrier and is the additional barrier for attachment of adatoms to the 
edge of a critical island (the limit where the island edge is a perfect sink is obtained by taking = 0). v is the 
attempt frequency, assumed to be the same for all microscopic processes. Thus S* = D exp(— /3E£), where D is the 
diffusion constant. 

N* can be expressed in terms of n via Eq. ([!]) , and thus the nucleation term, /, is a function of the average adatom 
density. To estimate n, consider the density of adatoms, n(r), around a typical stable island of radius R (r is the 
distance from the center of the island and radial symmetry is assumed). Under conditions of complete condensation, 
where no evaporation occurs, n(r) obeys the diffusion equation 
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d 2 n(r) 1 dn(r) 
dr 2 r dr 



F = , (5) 



where the quasi-static approximation has been used, suppressing the time derivative of the adatom density. This 
approximation is valid when diffusion is fast enough so that at each instance of time, n(r) reaches a quasi-steady 
state, where the flux is almost entirely compensated by the attachment of adatoms to stable islands. The relevant 
solution of Eq. (pf) obeys the following boundary conditions: 
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r=R - S 
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where S — D exp(— (3Eb), Eb is the additional energy barrier for attachment of adatoms to the edge of a stable island, 
and L is half the distance between stable islands (N s — 1/ttL 2 ). The first boundary condition holds for large islands, 
and was derived by Bales and Zangwill . 
The solution of these equations is 

FL 2 r F 2 2N FVtiL 2 D-S{^ R 2 \ 
Following Stowell and coworkers JL2| , it is easy to obtain an expression for n: 



n(L 2 - R 2 ) 



2nrn(r)dr 



R 



(8) 



F / i a 3\ 1 , F D-S / 



where it is assumed that the density of stable islands is already large enough so that « R 2 /L 2 and the contribution 
of n to the coverage is negligible. Since this theory is valid only in the small coverage limit, terms that vanish when 
9 — > were omitted. 

Eq. (Q) can now be rewritten as 

dB p ' W 

where the relation 9 = Frit was used, n is given by Eq. (||) and E* = Ei* . 

Eq. (|^) can be solved numerically starting from some initial condition N s (F,T,9i) = Ng 1 ' 1 (0i ^ since the quasi- 
static approximation does not hold when 9 = 0). But it is useful to consider two limiting cases where the problem 
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can be solved analytically. These limits occur when one of the two terms in expression for n is small and can be 
neglected. More quantitatively, the ratio of the two terms is 

G(F,T,0)= - 2 ^ (e 0Eb - l)y/N~ a . (10) 

If G is much smaller than 1 (limit I) for 9 > $i, the second term on the r.h.s. of Eq. (|§|) can be neglected. This is 
the perfect sink limit S — ► D (Et, — ► 0), and in this case standard results 0] are expected to hold. In limit II, on 
the other hand, G 3> 1 for 6 > $i, and the second term in (||) dominates; this is the case where exp((3Eb) ^S> 1, and 
the island-edge barrier is most important. In limit I, attachment of adatoms to island edges is infinitely fast and the 
kinetics is diffusion limited, whereas in limit II diffusion is fast, and the kinetics of adatoms is limited by the slow 
attachment to island edges. In both limits, Eq. (^J) can be written in the form 

Taking for simplicity E£ — Ej,, one gets in limit I 

T = i* + 2 , S=E*+i*E d 
and 



H{9) = 



lnfl-3/2^' +1 (12) 



while in limit II 



r= il+3 j £ = E* +i*{E d + E b 



~2~ 

and 



(13) 



The solution of Eq. (O) is 



Assuming <C N s (8), the second term on the l.h.s. of (|14| ) can be neglected and N s takes the form 

N s (F,T,6) = [Ta J\(e')d6f/ r e r 3£ / r /F . (15) 

Therefore, in limit I 

N s ~ e P(E'+i*E d )/(e+2) pi-/(i'+2) ^ ^ 16 ^ 

while in limit II 

N s - e 2 ^ E '+^( E d+E b )]/(t'+3) F 2t'/(t'+3) _ ^ 



(N s Q) r - {N^nf = Tae 0£ j" H(9')d9' . (14) 



The coverage dependence of N s also differs in the two limits. However, it is known |18|,|13|] that the methods that have 
been used in this work are not suitable for an accurate calculation of this dependence. For this reason, only the flux 
and temperature dependence of N s are emphasized in Eqs. ( |l6| ) and (p^). 

As expected (see above), in limit I the result (Eq. (|l6|)) coincides with the standard result calculated under the 
assumption that island edges are perfect sinks for adatoms. Eq. ( |T7| ) clearly shows that in limit II, when island-edge 
barriers are important, the behavior of N s as a function of F and T is strikingly different [ jlSj ] . The most unambiguous 
information about the importance of island-edge barriers in a specific experimental system can be obtained from the 
functional dependence of the island density on flux. It is a simple power law, iV s ~ F x , with an exponent X that 
depends only on the critical island size i* . While in limit I the exponent is in the range 1/3 < Xi < 1, in limit II it 
can be larger than 1 (1/2 < X\\ < 2). For a given value of i*, X\\ is significantly larger than Xj, a difference that 
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can be detected experimentally. If one has some information about i*, a measurement of N S (F) can indicate which 
of the two limits is more appropriate for the experimental system in question. Furthermore, if there is no a priori 
knowledge of i*, one can still identify the relevant limit (and thus evaluate i*) if X < 1/2 or X > 1. The former case 
is possible only in limit I, and the latter only in limit II. Once the proper limit has been identified and i* evaluated, 
the temperature dependence of N s can be used to estimate energy barriers; in particular, one can evaluate £. 

Which experimental systems are suitable for such a study? Homoepitaxial growth experiments of semiconduc- 
tors and metals with and without surfactants may be good candidates. For example, consider the experiments of 
Voigtlander and Zinner |2(J on submonolayer growth of Si/Si(lll) with Sb as a surfactant. In this case, the addi- 
tional island-edge barrier Ef, corresponds to E ex — Ed, where E ex is the barrier for exchange of an adatom with a 
surfactant atom at the edge of a stable island, and Ed is the barrier for diffusion on top of the surfactant layer. Limit 
II corresponds to the model proposed by Kandel and Kaxiras to describe surfactant mediated film growth; they 
assumed that surfactant atoms passivate island edges. Limit I is associated with the more standard approach to the 



same problem 21 , where no island-edge passivation is taken into account. An experimental value of X > 1 would 
therefore indicate that the model of Kandel and Kaxiras is adequate for this system. A value of X < 1/2, on the other 
hand, would favor the standard approach to surfactant mediated epitaxy. Intermediate results would be inconclusive. 

It is instructive to use experimental numbers in Eq. (10) for the function G, and estimate what Ef, should be in 
order for the system to be in limit II. For example, at T — 900-fT Voigtlander and Zinner measured an island density 
of N s pa 10 11 cmT 2 in the case of Sb mediated growth (see Fig. 3 of Ref. pofl ). The coverage was 0.15 bilayers, and 
since fi pa 10" 15 cm 2 , the estimate of G is G pa 0.23(exp(/3-E b ) - 1). Therefore, for E b = 0.4 eV G pa 39, and the 
system should be in limit II. Limit I would be clearly observed only if Eb < 0.1 eV. 

Another work of relevance here is the experiment of Andersohn et al. ]2^] , who measured the exponent X for 
homoepitaxial growth of silicon with molecular beam epitaxy (MBE) and CVD using disilane (S12H6). They concluded 
that in MBE X ps 0.75, which corresponds to i* = 5 — 7 assuming that the standard limit I applies. In CVD, on the 
other hand, they obtained X pa 1.25. As they emphasize, standard rate equation theory cannot explain this result, 
since it always predicts X < 1. In the framework of the present work, the result X = 1.25 is a natural consequence 
of a significant island-edge barrier (limit II) and corresponds to i* = 5. In fact, since during CVD with disilane, 
hydrogen is always present on the surface of the growing film, it is possible that hydrogen atoms bind to island edges 
and have to be displaced before a silicon atom can attach to the island. In this sense, hydrogen acts as a surfactant 
in this system pd| , and the present interpretation of the experimental result favors the mechanism of Kandel and 
Kaxiras over the standard surfactant mediated growth mechanism in this particular case. 

Here again it is possible to use the experimental numbers in Eq. (|l0|). The CVD experiments at T = 800K show 
that for different values of the flux, the island density at coverage of 0.15 bilayers changes between 10 11 and 10 12 cm~ 2 
(see Fig. 4 of Ref. p3])- Taking the smaller of the two values with E b = 0.4 eV, one gets G pa 75, and the system 
should exhibit a limit II behavior. Thus, the analyses of both experiments shows that the values of the additional 
barrier, J3&, need not be rediculously large for the system to be in limit II. 

It should be emphasized that the CVD experiments interpreted above are not a reliable test of the theory due to 
the complexity of the processes involved in CVD. Hopefully, this work will stimulate experiments on simpler systems 
that will enrich our understanding of the role of surfactants in thin film growth. 

In summary, this paper shows that island-edge barriers have a profound effect on the density of 2D islands during 
submonolayer epitaxy, and as a result on the morphology of the growing surface. The novel behavior of the island 
density, predicted in this case, has been shown to be relevant for experimental systems, including surfactant mediated 
epitaxy and chemical vapor deposition. 
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